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providing information about the properties of the compactification of the extra dimensions. 
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1. Introduction 

By now, inflation |2|, ^ is well established by the observational cosmic microwave back- 
ground radiation (CMB) data |Q, |5|. However, the inflationary universe is more like a 
paradigm than a theory, since a specific compelling model is still missing. The best chance 
to come up with detailed working models is within a fundamental theory such as string 
theory. String theory realizations of the inflationary universe may be separated into two 
types, depending on whether the inflaton is a closed or an open string mode. A prime 
example of the former type is when the closed string mode is a modulus, while for the 
latter type it is when the inflaton is the position of a brane inside the compactified bulk 
in a brane world. This is sometimes referred to as brane inflation [^, and the simplest 
speciflc realization as the KKLMMT scenario 0. We would like to make a detailed com- 
parison of the predictions of this simple specific realization of brane inflation to the current 
cosmological data. This scenario offers potentially distinctive stringy signatures that may 
be detected. 
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Brane inflation proposes a string theory motivated mechanism for inflation that appears 
to be generic in weU-studied models which dynamically compactify the extra dimensions 
in a consistent way. It can have the additional nice feature of a natural ending when the 
branes collide, where the collision itself is useful for reheating and the possible production of 
cosmic strings. String theory dictates both the dynamics of the inflaton and its potential, 
so that one can make precise cosmological predictions from a given set of background 
parameters. Furthermore, given the potential of current and future measurements of the 
power spectrum including tensor mode perturbation, non-Gaussianities and possibly cosmic 
strings, the string theory parameters can eventually be over-constrained by the data 1^. 
The calculability and limited number of parameters make brane inflation an interesting 
arena to explore the possibilities for cosmology in string theory. 

In this paper, we make full use of current data to constrain the model. The particular 
brane inflationary model studied here is a simple realistic scenario with the inflaton mass 
as one of the free parameters. As our understanding of the compactiflcation improves, one 
can probe numerous variations as well as more detailed features of the model. For this 
analysis we compare the model to the likelihood space allowed by Wilkinson Microwave 
Anisotropy Probe (WMAP) 3 year CMB data ||, §, 0] in combination with the Sloan 



Digital Sky Survey luminous red galaxies (SDSS LRG) matter power spectrum and 



the Supernovae Legacy Survey (SNLS) supernovae sample [|13 |. 

Consistency checks will provide strong limits on input parameters, and we can further 
narrow the choice using the measurement of the spectral index n^, the tensor /scalar ratio r, 
and current limits on non-Gaussianity. However, this data still leaves us with a few distinct 
regions of parameter space that we would like to examine in more detail. We outline below 
the basic differences between usual slow-roll and the brane inflation scenario, our results for 
calculating simple cosmological observables, and the features we expect in different regions 
of the input parameter space. If brane inflation is correct, one may hope that data will 
rule out everywhere except for a point (or a tiny region) in the whole parameter space. 
Hopefully, this work also gives a clear idea how further theoretical analysis, together with 
more data, will likely over-constrain the model and provide a critical test of the brane 
inflationary scenario. Within the brane inflationary scenario, the cosmological data relates 
directly to properties of the compactiflcation of the extra dimensions. In contrast to today, 
where the standard model branes are sitting at a corner of the bulk, the mobile branes 
in the early universe probe a different part of, and complementary information about, the 
bulk space. That the inflaton is an open string mode and so takes the Dirac-Born-Infeld 
kinetic form in its effective action leads to interesting new features |14]. 

We shall present our analysis in three steps: 

1. We re-examine the analysis in the literature and flnd agreement of the model with 
data for some ranges of the inflaton mass m, which is treated as a free parameter 
here. There is (i) the small m region @, IH], which is essentially the KKLMMT 



model, (ii) the intermediate mass region, where the tensor mode can be very large 
and which deviates from the usual slow-roll relation between the size of the tensor 
mode perturbation and its spectral index nt |17], and (iii) the large mass region. 
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where non-Gaussianity (due to the DBI action) may be large and detectable [14, 18 1. 
We find there are regions with a red tilted power spectrum, in agreement with the 
data in each of these cases. 

We consider the impact of requiring self-consistency, that the throat should be inside 
the compactified manifold, which must be finite for the Newton's constant to be 
non-vanishing. Requiring the size of the throat to be smaller than the volume of 
the bulk imposes a very strong condition on the size of the throat and so on how 
far the D3-brane can be away from the bottom of the throat. This condition, which 
is absent in field theory models, rules out eternal inflation of the random walk type 
in brane inflation This bulk volume bound also limits the e- folds allowed. In 

generic situations, imposing this bound here rules out both large tensor and large 



non-Gaussianity in the model |20|. This leads to the conclusion that the original 
KKLMMT scenario (with massless inflaton, very small tensor and non-Gaussianity) 
seems to fit the data best. 

3. We then discuss how the bulk volume bound may be satisfied in a number of variations 
of the simplest brane inflation model. As specific possibilities, we discuss a couple of 
scenarios where the large tensor mode may be present. Since this tensor mode would 
have a spectrum that deviates from that in the slow-roll scenario, this may be a way 
to pick up a string theory signature. These scenarios also allow large non-Gaussianity. 

The paper proceeds as follows: in section ^we outline the brane inflation model. In sec- 
tion ^ the background and perturbation evolution and the dependency of power spectrum 
observables on brane/throat parameters are discussed. Four different inflationary regimes 
in the DBI model are discussed. Section ^ presents the constraints on the inflationary 
model in light of cosmological observations. In section |5|, we discuss a specific scenario 
where the tensor mode can be large. In section ^ we tie together our findings and discuss 
implications for the future. Unless otherwise specified, mass units are given in term of the 
reduced Planck mass. 



2. The brane inflation model 

Suppose our universe today is described by a brane world solution in Type IIB string theory. 
In this scenario, 6 of the 9 spatial dimensions are dynamically compactified to a finite size. 
There are Z)3-branes (and probably D5- and D7-branes as well) that span our observable 
universe. Standard model particles are the lightest open string modes which must end on 
branes. The D3-branes are point-like in the 6-dimensional compactified manifold known 
as the bulk. Closed strings can be everywhere, with the graviton being the lightest mode. 
This fiux compactification that is dynamically stabilized automatically introduces a warped 



geometry, with regions in this bulk known as warped (deformed) throats |21, |2^ . 

In the early universe, it is easy to imagine additional Z)3-i)3-brane pairs around. A 
Z)3-brane has the same tension but opposite Ramond-Ramond (R-R) charge as a D3- 
brane, and is attracted (strongly) to the bottom of a throat. Inflation takes place while 
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the mobile D3-brane is moving towards the Z)3-brane and inflation ends when they collide 



and annihilate each other 24, 25 1. Fluctuations that are present before inflation, such 
as defects, radiation or matter, will be inflated away. Here, the Z)3-brane position (j) is the 
inflaton and the inflaton potential V^(0) comes from the brane tensions and interactions. 
The I?3-l)3-brane annihilation releases the brane tension energy that heats up the universe 
to start the hot Big Bang epoch. So brane inflation relies on the dynamics of the L'3-brane 
extended in our usual space-time and moving in the extra dimensions. Together, the 10- 
dimensional metric takes the form: 

ds^ = h^{r){-dt^ + a{tfd^^) + h-'^{r){dr^ + r'^dT.%J, (2.1) 

Here the cosmic scale factor a{t) is that of an expanding homogeneous isotropic universe 
spanned by the 3-dimensions x, and r is the coordinate along the throat. The warp factor 
h{r) and the metric on the extra dimensions (dS^^) are inputs to the model. Warped 
spaces are natural in string theory models and are useful for flattening potentials and for 
generating a hierarchy of scales with the UV at the top (edge) of the throat and the IR scale 
at the warped bottom (around r ~ r^). We will use one particularly well-motivated form 
for the full metric with h{r) ^ r/R, where i? ^ r^i is the scale of the throat. The throat will 
be described by these two numbers giving its length and warping. The inflaton </> is related 
to the position of a space-time filling L'3-brane moving in such a throat. Specifically, 



Tsr (2.2) 

where = [{2it)'^ gsa'"^]^^ is the tension of a D3-brane. Here Qs is the string coupling and 
the Regge slope a' = m"^ sets the string scale, where is the string mass scale. 

String theory suggests not only the metric for this model, but also the appropriate 
action and potential. Besides a topological piece (the Chern-Simons term), it is well- 
known that the world volume action for a D-brane also involves the Dirac-Born-Infeld 



(DBI) action, which has an unusual kinetic term |26, 27 1. Here we will use a particularly 

(2.3) 



simple background geometry, so that the action is correspondingly simple |T 

d^xa^(t) T^l-(t)'^/T + V{cp)-T 
where 

T{(p) = T^h\(p) (2.4) 

is the warped Z)3-brane tension at (j) and V{(j)) is the infiaton potential. The DBI action 
has been used in string cosmology before p8| , |29[| . The important new ingredient here is 
the combination of warped geometry with the DBI action, a consequence of the dynamical 
moduli stabilization in flux compactification. To see the key point, it is natural to introduce 
the Lorentz factor, 7, that tracks the motion of the brane : 

7(</.) = , ^ (2.5) 
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Note that the inflaton speed is hmited by the warped tension, i.e., cf)^ < T{(f)) and r((/>) 
decreases towards the bottom of the throat, irrespective of the steepness of the inflaton 
potential. As the limiting speed is decreasing rapidly down the throat, this can lead to a 
scenario where the inflaton is moving slowly but ultra-relativistically, generating predictions 
that are quite different from those of the usual slow-roll case. In the non-relativistic limit, 
7^1, the model ( |2.3| ) reduces to those with a standard canonical kinetic term. 

While the form above is useful for the calculations that follow, it is worthwhile to 
display the more complete and fundamental expression for the action 1 30]. For interesting 
and straightforward generalizations of the background geometry there can be many fields 
appearing non-trivially in the action, including the dilaton the metric G^u^ the anti- 
symmetric tensor B^u and the gauge field F^jj. The Chern-Simons term contains couplings 
between the brane and R-R fields (p-forms) Cp, with p even for type IIB theory. We use 
variable ^ and indices {a, 6} for coordinates and quantities on the brane. Then, ignoring 
the potential which we will discuss in detail below, we have 

5d3 = -Ts j d^ee"* yJdet\Gab + Bab + 2^a'Fai\ (2.6) 



Ma 



p=0 



ATr 



where the quantities inside the determinant have been pulled back onto the brane and 
the lb is for a brane/anti-brane. In particular, note that the pulled back metric (Gab) will 
contain the warp factor. In the second term ^us = gsTs is the brane charge. The integration 
is over the D3-brane world volume, and contributions should, of course, have the correct 
dimension. In the simplest cases, many of the fields appearing in ( |2.6D are trivially zero. 
However, for more interesting solutions we will need to solve the supergravity equations to 
find each of the fields. 

We see now that the simple case in ( |2.3D has constant dilaton (set to 0), vanishing 
pullbacks Bab and Fab, and only C4 non-zero (D3-branes are charged under the 4- form 
RR field C4). The Z?3-brane alone in this background is supersymmetric. In addition, 
the supergravity equations require that components of C4 on the brane be h'^{r)/gg. We 
have aligned the brane coordinates with the usual space-time coordinates so that the only 
non-zero derivatives in the pullback are those with respect to time, and we assume that 
only the radial motion of the brane is important. Then (|2.6D simplifies to 



Sd3 = J d^x a\t)T^h\r) (^-^l - h~^P + l) (2.7) 

From the supergravity solution, we find that the throat should be characterized by the 
parameters Ha = h{rA)-, the warp factor at the bottom, of the throat, where r = va, and 
a background number of charges we label Na- The scale of the throat R is given by [^] 

i?^ = A7rgsNAa'^/v (2.8) 

where v is the ratio of the dimensionless volume of the space X5 (in (p.lj)) to that of a 
5-sphere (which has volume vr^, assuming unit radius). For the well studied Klebanov- 
Strassler (KS) throat, which has a simple geometry (to be discussed in section |3.3| ), v = 
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16/27 [^l|. However, the KS throat has an S"^ symmetry, so one can consider its orbifolded 
versions, where v can be quite small (<C 1). We sketch the throat in Fig.||. 




Figure 1: A cartoon of the throat in the compact extra dimensions. There may in general be 
warped throats of a variety of sizes attached to the bulk space. R sets the scale of the throat, while 
hA is the warping at the bottom. The _D3-brane moves down the throat, attracted by a Z)3-brane 
(or a stack of them) sitting at the bottom. The inflation, 0, is related to the branc position, r. 



The inflaton potential can in principle be computed from a detailed knowledge of the 
string theory background (the flux compactification) and is of the form 

V{4>) = ^mV' + Vo + Vj,Q + ... (2.9) 

Let us first consider the constant and brane/anti-brane interaction terms, which in more 
detail are (for 7 ~ 1), 

Vb + Vj^r, = 2T,h\ f 1 - -L^l = ^0 f 1 - T^-^l (2.10) 



The constant term, providing the vacuum energy for inflation, comes from the warped 
tension of the brane and anti-brane. The interaction term comes from the attraction 
between the mobile Z)3-brane and a Z)3-brane sitting naturally at the bottom (small r 
end) of the throat. Its form can be understood in several ways, most simply as the analog 
of the Coulomb interaction in the 6 dimensions transverse to the branes. The interaction 
can be explicitly derived in flat space by examining closed string exchange between the 
branes, and the calculation can be extended for relativistic branes to determine 7 factor 
corrections (see ( |3.45 )). Finally, in the warped geometries considered here, one brane can 



be treated as a probe moving in the compact space with a metric perturbed by the other 
brane This explains the factor of 1/Na-, since the perturbation due to a single extra 
brane should be suppressed by the number of background branes. Note that the singularity 
of the Coulomb term is never reached. When the branes are within the string scale distance, 
the Coulomb potential is curtailed due to the exchange of the infinite tower of closed string 
modes, as well as the emergence of a tachyon mode, whose rolling signifies the collision and 
annihilation of the D3-D3 brane pair. 
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So we have 



(f>e>(f>i> 4>pivot > 4>E > 4'A 



where (/>e is the value at the edge of the throat, (f)i is at the onset of inflation, (ppiuot is the 
position when the observed CMB scale was originally produced, (/>£; is when the tachyon 
emerges and inflation ends, and <j)A at the bottom of the throat. 

While the constant and the coefficient of the Coulombic term are easily calculated in 
most cases, finding the effective mass, m, is quite involved [32, 33]. Contributions from 
the Kahler potential and superpotential after moduli stabilization may be large Q, and 
some work on simple models indicates that a' corrections are also important |34, 35, |36|| . 
In realistic examples, calculating m is a long and cumbersome task (although the situation 
is improved somewhat by the techniques of Because of these difficulties, we take the 

pragmatic approach of considering m to be a parameter of the model. The cosmological 
predictions are very sensitive to the value of and actually the shape of the throat 
pO| , [37[| . In fact, the size of m partially dictates how closely this model mimics standard 
slow-roll and where it deviates. Furthermore, as shown in Ref.[^], even a small in the 
slow-roll case modifies the power spectrum from red tilt to blue tilt. Since the value of m? 
is dictated by the details of the compactified manifold, CMB data could reveal to us the 
structures of the manifold, a very exciting aspect. There may also be other terms in the 
potential. For large values of (/>, we should include additional terms like (/>^, which we will 
discuss later in section 

Cosmic strings are generically produced during the brane annihilation towards the end 
of inflation. Although we shall not discuss the cosmological implications of cosmic strings 
here, we shall occasionally quote the prediction of the cosmic string tension /j, in terms of 
the dimensionless parameter G^. 



3. Cosmological implications of brane inflation 

In this section we outline the background expansion and the spectrum of density and tensor 
perturbations for the DBI inflation model. 

3.1 Homogeneous background evolution 

There is an important difference between the DBI scenario and the more usual slow-roll 
models often presented, which shows up in the kinetic term for the inflaton. The important 
feature of the DBI action (^j^) is the effective speed limit imposed by the Lorentz factor 7 
dU) introduced above [|8|, ||]. 

Notice that the speed limit depends on the warping, so that deep in the throat (small 
r and 0), 7 may grow quickly. This effect allows the model to achieve enough e- folds 
of inflation even in a steep potential and is the source of several interesting observable 
differences from the standard slow-roll. We emphasize that this behavior is an unavoidable 
feature, coming from the fact that the brane position is an open string mode. 
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Using the action in ( |2.3D modifies the usual cosmology [14|. The pressure p and density 
p are given by 

p= r(7-i) + y (3.1) 

p = r(7 - l)/7 - V 
The equation of motion for the inflaton is 

<A - + T' + ^0 + {V - T') ^ = (3.2) 

where 7, H, and V are all functions of (p and the prime denotes derivatives with respect to 
0. We have used the continuity equation 

P = -'iH{p + p). (3.3) 



Notice that for small velocity, (|3.1|) and (|3.2D reduce to the usual expressions. Some algebra 
gives a useful expression for (j) 

-2H' , , 

(3.4) 



We emphasize that we have not neglected the (\) term. From the Friedmann equations, we 
need w = p/ p < —1/3. Using (3.1), 

a V r(7 + 2 - 3/7) 



a 3M2 6M2 



(3.5) 



where the kinetic term always contributes negatively. Then 

2,MlH{ct>f = V{(p) + r(7((/>) - 1) (3.6) 
Using the solution to the equation of motion for 0, we may write 



-f{(l)) = ^l + m^T-^H\<P? (3.7) 
Physically, this model has a sound speed that changes with 7 as 

Cs = - (3.8) 

7 

The sound speed plays a role in the scalar power spectrum, discussed below, and in the 
shape of non-Gaussianities |39|. Equations ( p. 6] ) and ( p.7D give a differential equation for 
H{(j)). Once the Hubble parameter and 'y{4>) are known, one can find expressions for the 
power spectrum observables. We use the following parameters: 

(3-9) 



2M2 




7 


\ m) ) 


2M2 


(H"m 


7 


{ H{<P) J 


2M2 




7 


\H{ck) 7 


e 


2e - 277 + K 
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The far right hand sides and fj are for comparison with [40| and |39]. Note that inflation 
takes place for < e < 1, since 

- = H\l - e) (3.10) 
a 

All three parameters must stay less than one for the expressions for the observables to be 
valid. In the non-relativistic limit 7^1, they are simply related to the usual slow-roll 
parameters esR and ijsr: e esR, r] 7]sr - esR and k^O, where esR ^ ^{V /V)'^ and 
VSR = V"/V. 

3.2 Scalar and tensor perturbations 



The scalar density perturbation has been studied in Ref.|41, 18, 42]. Decomposing the 
inflaton into its rolling background (j){t) and a fluctuation 6, 

^ = (l){t) + S{x,t) (3.11) 

with a scalar perturbation (the Newtonian potential) $ in the de-Sitter metric, 

ds"^ = -(1 + 2^)dt^ + a{tf{l - 2^)dx^ (3.12) 

it is easy to see that 5 and $ are related so that there is only one independent scalar 
fluctuation. The linear combination 

becomes frozen as it exits the horizon during inflation, later generating the temperature 
fluctuation in the cosmic microwave background radiation. The evolution of C, obeys a 
linearized Einstein equation. In terms of the variable 

z = (3.13) 

one introduces the scalar density perturbation u = C,z. Introducing the conformal time r, 
dr = dt/a, we see that as a function of the wavenumber k satisfies 

+ 7^-73Zo Wfc=0 (3.14) 



Note that the fluctuations u travels at the sound speed Cg, 

cl = ^/^ = A (3-15) 

so Uk freezes when k crosses k = aH^ instead of /c = aH. To solve for Uk-, we express the 
potential in terms of inflationary properties. Following ( 3.13| ), we obtain 



= a^H^W (3.16) 

z dr^ 
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where we introduce the derivative of the DBI parameters with respect to the e-fold number 
N, e.g., rjN = V,N= 



H 



H 



where 



—En = e(2e — 27y + k) 

-VN = r/(e + k) - ^ 

—kn = k{2k + e — r]) — euj 

4Mp4 f H'{cP)H"'{cP) 
2M2 /y 



(3.17) 



(3.18) 



7 V 7 

Expressing the Hubble scale in conformal time (valid when e is roughly constant), aifr(l - 



1, ( 3.14 ) becomes 



1^ + I C'^^ 



1/4 



Uk = 



(3.19) 



where 



W 



+ 



1 



Here, 1^ 3/2 as the DBI parameters vanish. Since all DBI parameters, as well as H, vary 
much more slowly than a{t), we may take u to be close to constant, so the above equation 
behaves as a Bessel equation. We see that, for aH^ ^ k, the growing mode behaves as 



\uk\ 2 

so the spectral density is given by 



V]l\k) 



fc3 

2^ 



z 



27r|0| 



(3.20) 



(3.21) 



k=aH'y 



Then the scalar spectral index is given by 



n. 



1 = 



din k 

~ (1 + e + K)(-4e + 2r]- 2k) 



(3.22) 



which reduces to the usual slow-roll formula — 1 = —6esR + '^VSR in the limit 7^1. 
As the Z)3-brane moves down the throat, 7, n and e tend to increase. In the large 7 limit 
in the AdS throat, k — > — 2e while rj ^ 0. Then — 1 — > 0. 
The tensor mode spectral density, to first order, is given by 

2H^ 



(3.23) 



k=aH 
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and the corresponding tensor power index: 



_ din Vh -2e 

This is non-vanishing even in the ultra-relativistic case. The ratio of power in tensor modes 
versus scalar modes is 

r = 1^ (3.25) 

7 

To keep r < 0.5, we would like 7 to increase as e does. However, the non-Gaussianity 
bound constrains 7 < 31. This bound can be saturated under certain conditions in the 
intermediate regime. In the intermediate regime, it is possible to have e ~ 0.2 and 7 ~ 1, 
so that r exceeds the current bound. 

Note that these equations suggest a key way in which one may distinguish usual slow- 
roll inflation from the DBI scenario: the consistency relationship between the tensor power 
index and the tensor/scalar ratio is modified 

= --S (1^) (^•^'" 

This is in fact a common feature of non-standard inflation models with sound speed Cs less 



than 1 |41]. To first order, nt = — r/(8cs). Since 7 is greater than one (and grows with 
time), finding the magnitude of rif larger than r/8 would be evidence for a non-slow-roll 
scenario. It was pointed out in Ref . [^] that the non-Gaussianity f^i is proportional to 
7^. The present observational bound from non-Gaussianity {/nlI — 0.327^ < 300 yields 
7 < 31, where the limit on non-Gaussianities coming from DBI-type models is discussed 
in Ref. p% |43[|. The calculation of the full bispectrum for models of inflation with small 



sound speed (including DBI) is done in Ref. |39f| . The trispectrum is studied in Ref.|44|. 
We note that the four observables /nl, 1^, nt, and Ug can be related to the four parameters 
7, e, K and r]. 

3.3 Impact of throat warping: The Klebanov-Strassler Throat 

The power spectrum index is sensitive to the deformation of the throat once we are away 
from the slow-roll region, and as such the precise shape of the throat, i.e., the warped 



geometry and the deformation can be measured using observations |45]. For this reason 
we discuss the warped throat in some detail here. 

Consider two different approximate ways to incorporate the deformation of the warped 
conifold, namely the AdS cut-off and the mass gap (MG) cases for which the D3-brane 
tension is given by: 

TAds{(l)) = y , 4>>(pA 



X = T3R 



2 + 62^2 

A 

4 _ 27 TVa 



TMGi^) = (3.27) 



16 27r2 
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We expect (j)A — b, so away from the bottom of the throat, for large (p, we expect httle 
difference in physics. However, in Ref.[17|, we find — 1 is slightly positive (~ while 



in the second case, Ref.|46] finds that the spectrum has a red tilt. This difference appears 
in the DBI region. We expect a small (hopefully observable) effect of the warp factor in 
the large mass region. (The warp factor drops out in the slow-roll region.) This means 
and its running will tell us about the shape of the throat, an exciting prospect. Note that 
the location and size of the 4-cycles in the bulk (partially) measures the inflaton mass pO| . 
So we really can learn a lot about the compactification geometry from the CMB. 

Now we have to consider the geometry of the throat more carefully. Although there is 
a wide class of possible geometries, we only know the full metric of the Klebanov-Strassler 
(KS) throat. So let us consider this case. Consider the lOd-metric ( ^.1[) , 

ds'^ = h{rf{g^^d''xd''x) + h{r)-^dsl 

= h{rf{-df + a\t)d^'^) + h{ry\dr^ + r^dS^) 

4/3 

= hir)\-dt^ + a\t)dy?) + /i(r)-2(_l__dr2 + ...) (3.28) 

(here, e is the deformation parameter, not the inflationary parameter) where 

K[r) = (^i"M2r)-2r)V3 ^^^^^^ 
2^/3 sinhr 

Unless we want to study the iso-curvature density perturbations, we shall ignore the shape 
encoded in dT?. (For slow-roll, we expect this to be small, but they may be not in the DBI 
case.) The warp factor /i(t) is given by the following integral expression [47| 

h-\T) = 22/3 {g.Ma' f e-^/^ j(^) ^ (3 30) 



r(r)= r dx "" "^"^^ „^ ^ ( sinh(2x) - 2x ) (3.31) 
Jt sinh X 



where M is a parameter characterizing the background flux. Both K{t) and /(r) are 
well-behaved functions of r: 

K{t ^ 0) ^ (2/3)^3(1 ^ l! ^ ^ 00) ^ 2^/3 e-/3 

/(r^O) ^ ao + 0(r2), i[t ^ 00) ^ 2, .2"^'^ {t e~^^/^ (3.32) 

where ao ~ 0.71805. Here we take e^/'^ to have dimension of length. So we have 

h{r) = h{T) 

g2/3 

dr = —= dr (3.33) 

^/QK{T) ^ ' 

Integrating this equation gives the relation between the throat position, r, and r. 
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From the explicit expressions in terms of r, we can obtain h{r) or T[(j)). For large r, its 
relation to r in dsg is given by 

=^ W^'""^"^'' (3-35) 
In the limit e ^ 0, finite r implies r — > oo, so 

dsl dr"^ + r'^dY? (3.36) 
which is simply the conifold metric. At r = 0, we have Ha = h[rQ) = tq/R and 

h{rof = h{r = Of = co^'— = e-^-^/3^«^ (3.37) 

gsM a' 



where cq = l/(2^/'^y^), so 



1 fe^/^\ 



where R is given by ( p.^ ) and A^a = KM. 

The edge of the throat is glued smoothly to the bulk of the generalized Calabi-Yau 
(CY) manifold. The bulk of the CY manifold is defined to have no significant warping, so 
that the warp factor in the bulk is essentially of order unity. Suppose is the point where 
the throat is glued to the bulk. We may estimate Tc using ( 3.34| ): 



2/3 rrc ^ 2/3 i-Tc i 

-dT ^ R-r^^— -rjT-TNdr (3.39) 



^e-/=^«(i?-ro)e-2/3^1^!^ + l 



as long as Tc >> 1, which is true for reasonable parameter values. This expression agrees 
well with the numerical value for which h{Tc) = 1. 

The relation between Tc, e and Ig is important to impose the constraint that the throat 
must fit inside the bulk. The physical size of the throat is given by 

I = 6-^2 6^/3 r dr ~ v^/. (3.40) 

Jo K{t) 

1 /6 

We have assumed that I <Vq , where Ve is the volume of the compactification. 

Which approximate warp factor is closer to the actual solution depends on details of 
how the throat is cut-off. To illustrate this point, we show numerical examples of the 
various expressions for the warp factor in Figures |^ and ^. Note that when M > \/Na, the 
mass-gap formula is much more accurate than the AdS expression. In particular, the AdS 
warp factor indicates a throat that is much too short (the point where the throat connects 
to the bulk is when /i(tc) = 1). However, when M < \/Na, the mass-gap expression is 
actually too fiat and the AdS expression is more accurate. This is especially true for small 
Na, where the slope of the full solution changes considerably in one warped string length. 
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For completeness, we also plot the log-corrected solution (see [47|), although we do not use 
it in our analysis since it deviates strongly from the actual KS solution at small r. 

One can understand these results by examining the initial slopes of the various expres- 
sions, accounting for the change of coordinates between r and r in the small r limit. The 
important constant in the AdS and mass-gap expressions is rg/e^/'^ ~ (-v/iV^/M)^/^. If 
this constant is large, the slope of the mass-gap expression decreases and that warp factor 
flattens out, while the AdS expression is less negative. If the constant is small, the AdS 
case has a very large negative slope initially, while the mass-gap case has a stronger r 
dependence. 

3.4 E-folds and initial power spectrum observables 

The observable power spectrum parameters are measured at a range of wave numbers, k, 
which in turn correspond to a certain time range before the end of inflation. It is more 
convenient to translate this to the number of e-folds A'^e before the end of inflation. This 
number of e-folds is given by the integral 



H dt 



1 



2M| 



H' 



(3.41) 



In order to calculate the e-fold number correctly, it is crucial to determine the position of 
the brane when inflation ends, i.e. (j)E- From (|3.10 ), we know that inflation ends when 
e > 1. But in most cases, e remains less than one and inflation ends for a different reason. 
In the slow-roll scenario, rj grows much faster than e during inflation and it is actually 
r] = 1 that sets (pE- For the DBI scenario, e is suppressed by I/7 and almost never grow 
larger than 1 at the end of inflation. In fact, in the DBI scenario, inflation ends when the 
proper distance between DD pairs is the string scale Va'. 




Figure 2: Comparison of the various warp factor expressions. The solid (red) line is the full warped 
deformed conifold expression. The long dashed (green) line is the AdS warp factor. The dot-dashed 
(dark blue) is the log-corrected expression, and the short-dashed (light blue) is the mass-gap. For 
this plot Na — 10^, M = 4000. The vertical black line indicates te, where the tachyon develops 
and inflation ends. 
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Using the metric (|2.l|) , the proper distance between the branes can be calculated as 

dr 



ds 



h{r) 



For the AdS throat, with h{r) = r / R, we get 

For the mass-gap warp factor, with h{r) = y^r^ + r^/R, we find 



I sinh 



'a' 



R 



where (pQ = ^/T^rQ and the anti-brane now sits at </> = 0. 
3.5 Four inflationary scenarios 

The full expression for the inflaton potential is dependent upon 7, 



V 



m 



Vo (7 + 1)' 



(3.42) 
(3.43) 

(3.44) 



(3.45) 




-Log[h*'] 



Zoom In 




-Log[h^] 



Zoom In 





Figure 3: As in Figure ^ showing warp factors for (top panels) Na = 10^, M — 800 and (bottom 
panels) Na — 10^, M = 10. The right hand panels show a zoomed in region. The vertical black 
line indicates te, where the tachyon develops and inflation ends. 
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where the Coulomb term is important in the non-relativistic regime, 7 ~ 1, but does not 
play an important role when the inflaton is relativistic, 7^1. 

The inflaton mass is proportional to the Hubble scale so we can parameterize it as 
~ PH"^. Now (3 controls whether the mass term or the constant term Vq is dominant in 
the potential and thus determines whether inflation is slow-roll or relativistic. In the slow- 
roll case, we also have ~ Vq, so that f3 may be treated as a constant. For the relativistic 
case, the potential is dominated by the mass term (so we no longer have ~ Vq) but we 
can still parameterize m? in terms of the constant term Vq as m? ~ /SVq. 

Overall, there are four scenarios, three for a Z)3-brane moving down a throat and one 
for a D3-brane coming out of a throat : 

1. /3 ^ 1, 7 ~ 1, the slow-roll case, when ~ 0; the m? = case is studied in 
Ref.P] and the more general small case is studied in Ref.|15, IC], where it was 



found that (3 < 0.05; that is, the range 0.05 < /3 < 0.2 is ruled out. Likewise they 
find, for small /?, ~ 0.98 -|- /3, the tensor to scalar ratio log(r) ~ —8.8 -|- 60/3, and 
the cosmic string tension log(G/i) ~ —9.4 -|- 30/3. As we shall see in our analysis. 
Us — 0.97 when /3 ^ 0. To explain the difference here, one should note that when 
/3 = 0, — 1 ~ — 5/(3A'^g) Q, and in our analysis, the minimum number of e-folds is 



not fixed, but rather depends on the scale of inflation through (B.ll). So the allowed 
range of /3 is now more subtle as it depends on Vq which sets the inflation scale in 
the slow-roll case. 

2. /3 ~ 1, 7 ~ 1 during inflation, but 7 increases to a large value towards the end of 
inflation; this corresponds to some intermediate values of m?. This case offers the 
best hope of a large tensor mode that has a tensor power index nt different from that 
predicted by the slow-roll scenario [17|. This scenario is not easily studied analytically 



and is more suited to numerical integration as we discuss further scenario in section 
I 

We can see the origin of the large tensor modes by combining the differential form of 
the e-fold expression (|1|) with (U) to give the usual Lyth bound expression 

1 d(j) 2MpH' ^ 



(3.46) 

where r is a function of <j). For small r, one can get many e-folds for a relatively 
small displacement of (j). For relatively large r, it is difficult in usual slow-roll to get 
enough e-folds for an initial (pi that is bounded (i.e., less than Mp). To agree with 
data, we see that large r can be achieved during inflation only if r drops rapidly as 
(p decreases, and inflation ends. This is achieved in the DBI model when 7 increases 
rapidly as the D3-brane moves down the throat (r ~ 7~^, since the parameter e is 
nearly constant). This means that this intermediate scenario can have large r, but it 
is inevitably accompanied by large non-Gaussianity on smaller scales. 

3. /? 1, 7 is large, inflation is ultra-relativistic throughout. This region typically has 



large distinctive non-Gaussianity |18, 39| and has been further analyzed recently P5| 



We shall go into some details of this scenario in the next section. 
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4. For completeness, we would like to mention another interesting case, namely, when 
/? < 0, when the inflaton mass is tachyonic. The scenario becomes the multi-throat 
brane inflation scenario proposed by Chen p8| , 49 1. The Coulombic term Vq/) is 
negligible and inflation takes place as the D3-brane moves out of a throat (before 
it moves into a second throat). For small tachyonic mass, this is simply a slow-roll 
model. This DBI inflation can happen when inflaton mass takes a generic value, 
\m\ ~ H. The distance the inflaton travels through during inflation, Acj) « HR^^/Ts, 
is always sub-Planckian. This model may be realized in a multi-throat compactiflca- 
tion starting with a number of antibranes settled down at the ends of various throats. 
These antibranes are classically stable, but can annihilate against the fluxes quan- 
tum mechanically via tunneling |50| ], producing the D3-branes, which then leave the 
throat. This model predicts large non-Gaussianity with the same shape as in the UV 
model. The difference is the running, /jvl ~ 0.036(3'^N^, that is, here /nl decreases 
with k, while, in the UV model, Jnl increases with k. The power spectrum index un- 
dergoes an interesting phase transition at a critical e-fold, from red (n^ — 1 ~ —A/N(.) 
at small scales to blue [ug — 1 ~ ^/Nf.) at large scales [49, 51 1. If such a transition 
falls into the observable range of WMAP/PLANCK, it predicts a large running of n<j, 
i.e., a large negative dns/dlnk in some region, but is un-observably small otherwise. 

3.6 Ultra-relativistic inflation 

In the large 7 case, the inflaton mass m is relatively big and one may neglect the other two 
terms in the potential ( p.45 ). To compare with the discussion in Ref.|]l4|, let us choose 
the deformation to be given by a mass gap, 

2 + 62)2 54 262 „ 

— ~ 1 fA^ 

A 



ri(0) = T(0) 



A A A " 

where 6 measures the deformation and A ^ Na- The density perturbation 

-10 _ fV' _ X ^ 



(3.47) 



P{k) ~ 10 



367r2M4 1447r2 



+ 62 



2 / X 4 

m 



Mr, 



(3.48) 



and 



so 5h — N^/V^, yielding 



A ~ lO^'^ — ~ 10~ 



2M„ 



m 



,2 1 



^2 + 6' 
6| + 62 



(3.49) 



(3.50) 



In terms of the throat length R, the value of the inflaton at the edge of the throat is 

(J)^ = y/j^R ~ m^A^/^ (3.51) 



and 



6 = miRhA 



(3.52) 
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For the SUGRA approximation to be valid, we require 

<mi ^ — <A-V4~10-^ (3.53) 



2 

Combining (|3.50|) and (p.53D, we find 



ra 



Mp 



> 10-2 (3.54) 



Now we have, for non-Gaussianity, 



where fig is the cosmic string tension. For Gfis < 10^^, /tvl will exceed observational 
limits. So the presence of large non-Gaussianity implies cosmic strings should be absent. 
This can happen if there is no l)3-brane at the bottom of the throat. 

If a l)3-brane is at the bottom of the throat, the tachyon appears at (pt = hATns, which 
is typically smaller than b. Also, we would like > 10~^ to ensure efficient (re-)heating 
fs^ . This choice is consistent with the other conditions so far. 

So far, we have not yet demanded that the throat must be inside the bulk. For typical 
throats that we know of, the width of the throat is comparable to its length, i.e., ~ R. 
If that is the case, then we find that 

^ < ^ 10-12 (3.56) 



which is clearly incompatible with the condition ( 3.54| ). Seeing non-Gaussianity (and as- 



suming a D3 brane inflaton), therefore, implies that the throat is not a KS throat or any 
of the Y^''^ type. The in the KS throat has an symmetry, so one may choose to 
orbifold it. Generically, we expect the size of the edge of the throat is R'^/p for Zp. For a 
large enough p, we may relax the above constraint ( ^.56 ). This constraint translates to the 



bulk volume bound (4J) introduced in Ref. [l53| since the infiaton cf) is measured in units of 
nis and the bulk size (which determines Mp) can be compared to the known throat scale 
i?^ oc NAa'^. 

As we shall see, consistent with this discussion, when combined with data, imposing 
the bulk volume bound would rule out the region of the parameter space with large 7 (and 
large non-Gaussianity). 



4. Comparison with cosmological observations 

In the previous section we outlined how the initial spectrum of scalar and tensor pertur- 
bations arising from brane inflation can be characterized in terms of six parameters: the 
inflaton mass, m, background charge, Na_, warping, hA_ = h((f)A), the orbifold/ volume ratio 
parameter in ( |2.8D , v, string length, and the string coupling, gg. To compare with 

observations we consider the predicted spectrum normalization, As{kpivot) = PR{kpivot)i 
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scalar and tensor spectral indexes, ns{kpivot) and nt{kpi^ot), and tensor to scalar ratio, 
r = Ph{kpi^ot) / PR{kpivot) , at a near-horizon physical scale kpi^ot = 0.002/Mpc. 

For the following numerical analyses, unless otherwise stated, we fix the string coupling 
Qs = 0.1 and length a' = 1000. We investigate the constraints on the remaining four 
parameters {NA,v,hA,'m'^}, and the slow-roll, intermediate and ultra-relativistic regimes 
outlined in section |3.5| . We generate Monte Carlo selected brane inflation models, obtained 
using a flat prior on {log N a, log v, log , As (kpivot)} and flnding Ha to give the required 
value of As^kpi^ot)- As mentioned in section |3^ , although the slow-roll and ultra-relativistic 
regimes can be easily studied analytically, the intermediate regime cannot. In appendix ^ 
we outline the numerical approach to integrating the equation of motion and perturbation 
equations, and in appendix ^ we describe how the inflationary initial power spectrum is 
translated to the physical scales today. 

We compare the model parameters, for both the AdS and mass-gap geometries, to 
observational limits on the amplitude, spectral index, tensor-to-scalar ratio for the initial 
power spectra from the WMAP CMB data |5|, ^, 10, 11], SDSS LRG matter power spectrum 



datafl^ and SNLS supernovae fl^]. We use the publicly available cosmomc code ||54| 
to perform a Monte Carlo Markov Chain analysis for a ACDM cosmological scenario, 
assuming flat priors on the following set of cosmological parameters: the physical baryon 
and cold dark matter densities, Uf, = Q^h'^ and = ^ch'^ (where h is the reduced Hubble 
constant h = Hq/WO), the ratio of the sound horizon to the angular diameter distance 
at decoupling, Og, and the optical depth to reionization, r, as well as power spectrum 
parameters, ns{kpij;ot), log As{kpi^ot)i and r(kpi^ot)- As we discuss below, and in flgures ^ 
and ^, in both the AdS and mass-gap scenarios, the running in the scalar spectral index 
is extremely small over observable scales, |(i?7-s/(iln /cl < 4 x 10^^, and the vast majority 
models closely obey the consistency relation nt = —r/8. We compare the Monte Carlo DBI 
models to compatible observational constraints obtained with dus/dlnk = and both with 
the consistency relation, and with r and nt allowed to be unrelated with —0.02 < rit < 0- 
We consider purely adiabatic initial conditions and impose flat spatial curvature. The 
MCMC convergence diagnostics is done on seven chains though the Gelman and Rubin 
"variance of chain mean" /"mean of chain variances" R statistic for each parameter. 
We also impose the following two consistency constraints: 

1. Sufficient e- folds and all e- folds in the throat: Although there may be inflation from 
other sources (other branes or other parts of the Wi trajectory) we require that 
sufficient e- folds to solve the usual horizon problem take place in a single throat. 
One may in principle loosen this requirement, but a more detailed knowledge of 
the geometry in the bulk space would be needed. The required number of e- folds 
depends on the reheating temperature Trh and inflation scale Hj, as outlined ( [B.ll| ) 
in appendix ^, roughly as ~ 68.6 -|- ln{Hi^ /Tp;j{)- This equation comes simply 
from the requirement that Ha at the beginning of inflation is at least as small as Ha 
today, where H is the Hubble parameter and a is the universe's expansion factor. In 
other words, the largest scale we observe today {kn ~ fepj^oi/lO) was once in causal 
contact. The expression is approximate because it assumes H is constant during 
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inflation. In the appendix we discuss how to implement the condition more precisely. 
The requirement that enough inflation occurs in the throat can then be expressed as 
(t>H < 4>UVi where (j)uy is the scale where the throat is glued to the bulk. We use 
the condition h{(j)ijv) = 1- For simplicity we assume 100% efficient reheating, so that 
Trh is set by the warped brane tension. 

A proportionally small throat: In string theory, the reduced Planck mass Mp is related 
to the (warped) volume VJj of the 6-dimensional bulk via Mp = V^/ {2'K)^T:g1a^ So 
the position of the L)3-brane cannot be physically larger than the largest dimension of 
the compactified bulk. For consistency of the string theory model, the throat region 
should be small compared to the bulk of the compactified space. Other aspects of the 
model (like moduli stabilization details or physics in other throats), therefore, should 
be insensitive to the details of the inflationary throat. Assuming that the throat has 
the geometry of the KS throat or its orbifold, one can find the bound on how far 
the inflaton can be away from the bottom of the throat. This turns out to be a very 



tight constraint, as shown by Baumann and McAllister |40|. One may put a minimal 
bound relating Mp to the parameters rus and Na by calculating the warped throat 
volume. The condition can be usefully expressed as a bulk volume bound on the field 
range of the inflaton , 

^p^^ot-<i^A (4.1) 



Mp V Na 

where (t)pivot is the value of cf) corresponding to the physical pivot scale today, kpiyot, 
where the normalization to CMB data is applied. 



Together with ( 3.46 ) and even considering that r may decrease rapidly during in- 
flation, Ref. |]40[| concludes that in DEI inflation with a quadratic potential, the ten- 
sor/scalar ratio would be unobservable for Na > 10^ and would badly violate the 
current bound on non-Gaussianity for Na ^ 40. One may rewrite the Lyth bound 
( |3l^ ) as 

f^fond / ^ \ 1/2 



V^CMB 

Note that this definition has number of e-folds M and Mes increasing with time. From 

(El), 

=^ = 2(. + «-,) (4,3) 

Then in the ultra-relativistic case, AAeflf ~ C'(e^^). In the most optimistic observable 
case and with Mes ~ 1, ^cmb ^ 10~^ requires Na < 10^. One must then check 
whether the current bound on non-Gaussianity is violated. For our purpose, this 



constraint in the form of the bulk volume bound (4J) turns out to be most useful. 
Note that this bound may be relaxed in some variations of the model that we discuss 
later. 
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Figure 4: Comparison of Ug in different warped geometries. 



We impose the first condition directly in the code, but implement the last as a consis- 
tency cut. Since other geometries (e.g., a squashed throat) or other choices for the inflaton 
may have a different dependence on A^^, it is interesting to examine the range of possible 
observables before and after imposing this constraint. 

In Figure Q we demonstrate the sensitivity of ns{kpiyot) to different warped geometries, 
holding all parameters, except P, at typical values and varying /? from slow-roll, /? ^ 1, 
to ultra-relativistic, /? ^ 1, regimes. We compare the AdS warp factor, the mass-gap 
warp factor and the exact Klebanov-Strassler warp factor. We see that in the slow-roll 
regime, since the warp factor does not come into play, different warped geometries give the 
same prediction on n^. When /? > 1 the model starts to deviate from slow-roll and the 
warp factor comes into play in the equation of motion, yielding different predictions on Ug. 
Typically, we find that the mass-gap geometry gives a more red-tilted spectrum than the 
AdS geometry and the exact KS is in between the two. 

In Figures |5| and ^ we demonstrate the wide variety of observational properties arising 
in DBI inflation. In addition to slow-roll behavior with red tilted, < 1, and blue-tilted. 
Us > 1, spectra, DBI inflation can also give rise to relativistic behavior 7 > 1, and interme- 
diate models with large tensors, small running and a strong running-spectral tilt relation. 
The figures show key 2D parameter spaces for both observables and model parameters 
with Monte Carlo generated models which are consistent with the WMAP+SDSS+SNla 
normalization (|B.1| ) at the 95% confidence level (c.L). The red squares are within la in the 
upper-left r — Ug panel, the orange-yellow triangles are within 2fT, while the black crosses 
are outside 2(t. The nine panels from bottom left are: 
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1. lo giQm? — logio^ panel. Low masses, < 10 ^'^Mp, produce slow-roll behavior 
7 ~ 1 with small tensors, r < 0.05, and obeying the usual consistency relation 
r = —Snt- Imposing a flat prior on logj^o"*-^ and log^^g'^) majority of slow-roll 
models have a red-tilted spectrum. The normalization and e-folding constraints allow 
the orbifolding parameter, v, to extend over a large range of values, ~ 10^^ orders 
of magnitude. The intermediate and relativistic regimes are constrained in a small 
section of this space with very small orbifold factors v < 10~^^. As such, only the 
slow-roll regime is consistent with the Klebanov-Strassler Throat for which v = 16/27. 

2. logio ~ logio panel. Two distinct regions, the slow-roll regime for < 
lO^^^Mp and a relativistic region with > 10- ^^M^ are carved out. The spread 
of Ha at small indicates the insensitivity of data to the warp factor of the throat. 
The requirement that there are sufficient e-foldings to include the physical scale kpi^ot 
gives the lower bound on Ha for a given mass in the slow roll region, with the blue 
tilted slow-roll models lying close to this boundary. In the intermediate and rela- 
tivistic regimes the normalization constraint translates into a thin region of allowed 
models. The relativistic models require a small warp factor at the bottom of the 
throat, hA ~ lO'^-^. 

3. log 

loi't'pivot ~ ^Pa) ~ log 10 -^A panel. Conveniently demonstrates the bulk volume 
bound ( [4.11 ). The relativistic models all have {(j)pivot ~ 

(t)A)/Mp > ^/4/Na, as such, 
imposing the bound rules out all but the slow-roll regime. 

4. dus/dlnk — Us panel. For both the AdS and mass-gap geometries the bulk of the 
models satisfying the e-fold and normalization bounds have a very small negative or 
near zero running in the spectral tilt, —0.001 < dug / dlnk{kpivot) < 0. The relativistic 
models, in particular, have extremely small running and ~ 1 as discussed in regards 
to (l3^ ). 



lo 

Sio fn^ — ns panel. For = 0.1 and a' = 1000 the slow-roll regime tends towards 
a common spectral tilt ~ 0.972 as m? — > 0, which reproduces the spectral index 
in the KKLMMT scenario 0. 

lo 

Sio 'm'^ — r panel. Models with significant tensor contributions with r > 0.05 are 
purely generated in the intermediate regime with m? ~ lO^^^Mp. 

Us — r panel. Models with a noticeable deviation from the Harrison-Zeldovich spec- 
trum (hs = 1 and r = 0) are possible with 0.96 ^ ng < 1.05, including models with 
significant tensor contributions < r < 0.15. ( 3.22| ) and ( 3.25 ) lead to the relation 



4(1 — ns)/7 ~ r. The AdS Monte Carlo models are able to have slightly higher tensors 
than the mass-gap geometry. The mass-gap, however, allows an interesting branch of 
models in the intermediate regime with large tensors and distinct vs. r relation. 

rit — r panel. Slow-roll and many intermediate tensor mode models are in good 
agreement with the r = —8nt consistency relation. Significant deviations from this 
relation are seen in the relativistic models where the models have a significant tensor 
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spectral tilt, but comparatively much smaller tensor amplitude given by r ~ nt/8j 
in ( p6D . 



9. 7 — r panel. There is a the clear division of behavior between the large mass models, 
with intermediate models possessing large tensor modes but a comparatively slowly 
rolling inflaton 7 ~ 1, and relativistic scenarios with 7 > 1 but vanishing tensor 
contributions in the limit 7 S> 1. For large 7, the Monte Carlo models are subject to 
additional constraints from CMB non-Gaussianity measurements, 7 < 31. 

In Figures ^ and ^ we compare the observational constraints from 3-year WMAP CMB 
+ SDSS +HST GOODs data with the Monte Carlo results for DEI inflation for the AdS 
and mass-gap warp factor respectively without imposing a constraint on (ppivot- The panels 
show the regions within la (dark blue) and 2a (pale blue) allowed by the data. The DBI 
models are marked red (yellow/orange) if they lie within la {2a) of the data in the key 
Ug — r space. Black points are outside the 2a region. It is clear that the data is constraining 
the parameter space of the model in a meaningful way. 

One can see that in both cases the inflationary parameter space {Na,v, /ia, m^} yields 
models wholly consistent with observations. For = 0.1 and a'Mp = 1000, the slow-roll, 
low mass, regime tends towards a common spectral index ns{kpi^ot) ^ 0.972, consistent 
with observations at the la level. An intermediate regime, for masses ~ lO^^^M^ and 
Na ~ 1, with a significant tensor contribution, is also consistent at la level. The relativistic 
(7^1) regime is consistent with data at the 95% c.l. For the AdS warp factor a handfull 
of Monte Carlo models, in agreement with the CMB and galaxy spectra at the 95% level, 
exceed the non-Gaussianity constraint from WMAP CMB on relativistic models, 7 < 31. 

In Figures |and|l^ we show the allowed Monte Carlo models after additionally requir- 
ing that they also satisfy the bulk volume bound ( |4.1[ ). This bound effectively constrains 



the tensor to scalar ratio, r{kpivot) < 4 x 10 In this case the slow-roll regime is allowed 



we 



but the intermediate and relativistic regime do not satisfy the cut. In Figures 11 
investigate the intermediate regime more fully, for a variety of specific parameter values, 
and find that the bulk volume constraint disallows any of the observationally consistent 
models to survive. 

For the intermediate region, it is worth noting that the initial (p takes a rather large 
value. In this case, we expect that the truncation, keeping only the m?'(j)^ term in the 
inflaton potential, may not be a good approximation. Let us, therefore, consider the 
impact of introducing a (f)^ term to the inflaton potential ( p. 91) : 

V{<t^) ^ y(0) + ^Fo i4A (4-4) 



We see in Figure |l^ that turning on a 0^ term tends to blue (red) tilt the power spectrum 
for small (intermediate) values of . 

In Figure |l5|, we show how the parameters and observables change with different values 
of /X. We find that as increases, Ka has to decrease in order to fit COBE normalization, 
meanwhile, Ug is driven to the red end and r increases. We also show in the figure that 
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including the (p'^ term does not help to reduce the value of (ppivot but makes it even larger, 
consistent with the bound ( ^.46 ) as r increases with larger //. 



5. Possible scenarios with a large tensor perturbative mode 

The intermediate region is severely constrained by the limited size of the bulk. If not for 



this bulk volume bound (4.1), we see that the tensor mode can be substantial. However, 
we would like to point out that there are ways in which future realistic models may in fact 
explore a wider parameter range. Let us discuss a couple of possibilities where the bulk 
volume bound may be relaxed enough to allow a large tensor mode. 

To see how to bypass/relax the bulk volume bound, it is useful to first recap the 
reason why it is so constraining. In a typical chaotic inflationary scenario, the inflaton 
field value can be much larger than the Planck mass Mp and the 2 slow-roll parameters are 
comparable, esR ~ rjsR- Since r = IQesR, we see why the tensor mode can be large (i.e., 
r is a few percent or larger). In brane inflation, the inflaton, which is the brane position, 
cannot be bigger than the bulk size, that is < Mp. To get enough inflation, it must move 



slowly. The Lyth bound ( 3.46 ) then tells us that r must be small if we want to get enough 
e- folds from moving down a throat. This also explains why esR <^ rjsR in this small field 
model. 

One way to satisfy the bulk volume bound is if we need only a few e- folds coming from 
the motion of (j) down the throat; a substantial number of e-folds are coming after cp has 
moved down the throat. This leads us to consider the following scenario. 

• Multi-brane-multi-throat scenario Perhaps the simplest way to generalize the single- 
throat scenario is to loosen the requirement that all inflation takes place when a single 
brane moves down a single throat. In general, there are a number of throats in the bulk in 
the flux compactification, with different warped factors. As suggested in ||5^, brane-flux 
annihilation at the bottom of one throat may lead to a number of n ^ Na D3 branes 
which can move out of that throat to the bulk. A number of these mobile Z)3-branes 
move in the bulk and then drop into these throats. So (j) is describing a specific Z)3-brane 
moving down a specific throat, say, the ^-throat. After it reaches the bottom, other D3- 
branes are still around and inflation continues. They can generate e-folds before dropping 
into some other throats, for example in some slow-roll scenarios. This implies that the 
CMB observables in WMAP correspond to a much smaller number e-folds away from its 
annihilation with the Z)3-brane at the bottom of the j4-throat. As an illustration, suppose 
there are 50 e-folds coming from other branes, then (j) has to go through only 5 e-folds before 



its annihilation. Following the Lyth bound ( 3.46| ), we see that typical r can be 100 times 



bigger than before, thus allowing a large tensor mode. This also requires a much smaller 
(f) to begin with, allowing the model to satisfy the bulk volume bound. The predictions at 
kn may then depend on more parameters. 

In fact, inflation can take place when the branes are coming out of a throat, in the 



case when m < 0, as in the Chen scenario [48, 49 1. This scenario will easily satisfy the 



bound ( ^4.11 ). However, a large tensor mode seems very unlikely since 7 is decreasing as the 
brane is coming out of the throat. 
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• D5-, D7- and/ or D9-hranes We have considered only the case of a Z)3-brane moving 
in the throat. While this may be the simplest and most natural scenario, one may also 
consider wrapped -D5- or Z)7-branes. Here, Z)5-branes can wrap a 2-cycle of the throat, 
while Z)7-branes can wrap any of the 4-cycles that are generically present in the flux 
compactification. A D9-brane wraps both cycles. In these cases, there are additional 
parameters involved coming from the wrapped cycles. The predictions of a wrapped Df>- 
brane depend on a different power of Na than the D2> case. It is easy to see that the change 
in predictions is related to the different normalization relating the canonical inflaton to the 
Db position [^5|. The Db action is proportional to the 1)5 brane tension, so once the 
integration over the wrapped two cycle (S2) is performed, the 4-dimensional action is 
proportional to T^Vols^- For a particular choice of embedding, one finds (A(/)/Mp)^ ~ 
y/gs/NAV. Here the orbifolding does not cancel out and a smaller power of Na appears. 
In addition, one may add fluxes on the brane which allow the model to match the data 
with relatively small Na- 

There are some other scenarios which, each by itself, will not generate a large tensor 
mode, but may be combined with the above scenarios and/or with each other to satisfy 
the bulk volume bound and enhance the magnitude of the tensor mode. 

• Squashed throat Although the warped geometries that are simplest to construct have a 
symmetry relating the size of the radial direction to the size of the cross-section, there is no 
reason in principle not to consider more general possibilities. A "squashed" throat would 
provide a different relationship between the inflaton range and the background charge, 
changing how the consistency bound is applied. 

• Motion at the bottom of the throat Recall that the throat is a warped deformed 
conifold, where the bottom of the throat has the geometry of a finite and a shrinking 
S'^. For simplicity, the scenario we have discussed so far puts both the l)3-brane at the 
bottom of the throat (at (^a) and the mobile D3-brane at the same position in the S^. In 
this case, when they are within a string scale distance (when at where (pE > 4' a), the 
open string tachyon mode appears, thus inflation ends before the Z)3-brane has reached the 
bottom of the throat. If the 2 branes are separated in the finite S"^, then one can get some 
number of e-folds after the Z)3-brane has reached the bottom but before it annihilates with 
the L'3-brane. This reduces the number of e-folds (p has to generate while it is moving 
down the throat, thus easing the constraint coming from the bulk volume bound. 

6. Conclusions 

We have discussed the cosmological predictions for DBI inflation, including the dependence 
of observed initial power spectra observables on the brane evolution and throat geometry. 
Three inflationary regimes are of interest, namely the slow-roll, intermediate relativistic 
and ultra-relativistic regimes, with rich properties including large tensor modes and non- 
Gaussianity. 

After numerically integrating the background and perturbation equations we compare 
the predicted power spectrum properties to current cosmological CMB, galaxy and su- 
pernovae data. We find that the three regimes can fit the data well although the power 
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spectrum normalization and e-fold constraints severely cut into the full multi-dimensional 
parameter space. 

The small rn? parameter region, which corresponds to a specific model within the usual 
slow-roll scenario, satisfies the cosmological data quite well predicting a spectral tilt around 
nsikpivot) — 0.972. Its predictions are in line with a typical infiationary scenario. 

Before imposing the bulk volume bound ( [4.1| ), we find that the (relatively) large 
region, which corresponds to the DBI ultra-relativistic scenario, is compatible with the 
present data, in agreement with previous analyses. The potential of having a large non- 
Gaussianity |18, ^] with distinctive bi-spectrum [^] and trispectrum [44| is very encour- 
aging to testing this very stringy scenario. However, imposing the bulk volume bound rules 
out most if not all of the large 7 region. That the cosmological data strongly restricts the 
allowed parameter space and brane inflationary scenarios means that the scenario can be 
tested; this is very encouraging. 

Interestingly, we find that the intermediate m? region, for which a numerical analysis 
is required, where the tensor mode can be large and the nt — r relation deviates from 
the usual consistency relation, is also consistent with the data, when one neglects the bulk 
volume condition ([4.1|). When this condition is imposed, however, the region is ruled out 
in the simple scenario. There are a number of ways to vary the model to relax this bound 
or the constraint coming from this bound; however, in spite of this, we do not expect the 
ratio r to be much bigger than a few percent. Two possibilities discussed here are the 
multi-brane multi-throat scenario and the wrapped D5-brane scenario. If r is large enough 
to be measured, the deviation from the r — nt relation in the slow-roll case will be a fine 
signature for stringy physics. 

The goal here is not only to find a stringy inflationary scenario that fits the data, but 
that such a scenario will also predict distinctive stringy signatures that may be detected 
and measured. Besides cosmic strings, searches for distinct string theory signatures in the 
slow-roll KKLMMT scenario continues. Recently, the proposal that the Seiberg duality 
cascade in the KS throat may lead to observable signatures in the CMB power spectrum 
1 53] is encouraging. Features in the warp factor or effects due to the fact that the infiaton 
is really a six-component field may source non-Gaussianity and slightly alter predictions 
for the power spectrum |5^, 57 1. 

Although the predicted tensor contribution of a few percent is beyond the capabilities 
of the WMAP satellite, this could well be in the detectable realm for the next generation 
of CMB instruments focusing on precision measurements of CMB BB polarization modes. 
Examples are, amongst others, the PLANCK, Clover, Spider, and BICEP projects, which 
hope to get Ar ~ 0.01 — 0.05. As such, the near future holds exciting prospects for better 
testing the rich properties of brane infiation. 

Finally, it is worth mentioning that cosmic strings are generically produced toward 
the end of brane infiation. We find that the cosmic string tension ^ roughly satisfies 
lO"^"^ < < 10^^; its particular value depends on the choice of m and A. Roughly 
speaking, small allowed values of G^, 1Q~^'^ < G^ < 10"^'' are accompanied by a red tilt 
in the power spectrum, and in some cases by large non-Gaussianity. 
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Figure 5: Taxonomy of the inflationary parameter space for the AdS warp geometry showing 
DBI inflationary models from a Monte Carlo simulation which satisfy the WMAP-|-SDSS-|-SNla 



normalization constraint [BA) at 95% c.l.. The flgure shows the wide variety of inflationary behavior 
arising from DBI inflation, including relativistic models (1.1 < 7 < 10,fllled yellow triangles and 
7 > 10 full red squares), large tensor modes (open green triangles), and blue and red-tilted slow- roll 
spectra (black and blue crosses, respectively). The taxonomy is presented in terms of relationships 
between the predicted spectrum observables (n^, r, dng/dlnk) and key model parameters {m^, Ha, 
Na, 7, and 4>pivot)- Note that the bottom right figure shows models in comparison to the bulk 
volume bound imposed by If the bound is imposed only slow- roll (low tensor, small running, 
non-rclativistic models) with (j>pivot — <t>A < v^T^^Ya are allowed. See the main text for further 
discussion. 
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Figure 6: As in Figure ^ but for the mass-gap warp geometry. 
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Figure 7: Inflation parameter space for tlie AdS geometry in comparison to cosmological observa- 
tional constraints from 3-year WMAP data in combination with SDSS LRG power spectrum data 
and SNLS SNla data. The top left figure (rig vs r constraints) and top center figure (n* vs r con- 
straints) of each set show the observational constraints without imposing the nt = — r/8 consistency 
relation, allowing —0.02 < nt < 0, at the 68% (dark blue) and 95% (light blue) confidence level. 
Overlaid are DBI inflationary models from a Monte Carlo simulation which satisfy the WMAP 
normalization constraint (BJ^) at 95% c.l. and are in agreement with the Us — r WMAP-I-SDSS 
LRG-hSNla constraint within 68% (red square), 95% (yellow triangle) and outside 95% (black 
cross) confidence limits. The constraints with the consistency relation imposed (dashed contours in 
the top left plot) yield a very similar splitting of the DBI models by confidence level. 
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Figure 8: As in Figure ^ but for the mass-gap warp geometry. 
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Figure 9: As in Figure [t] but imposing the bulk volume bound Q of (jipivot — (l>A < ^J^I^A for 
the AdS warp factor. Note that the tensor ratio and tensor tilt scales are reduced in comparison 
to Figure 0. The rescaled plots show that the, much more restricted, set of models are all in the 
slow-roll regime with r < 0.003 and 7 — 1 < 10^^ and satisfy the rit ~ —r/& consistency relation. 
The intermediate regime, while in strong agreement with observations, is incompatible with the 
bound. 
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Figure 10: As in Figure y but imposing the bulk volume bound for the mass-gap geometry. Only 
slow-roll models remain consistent with the data at the 95% confidence level with r < 0.001 and 
7 — f < fO^^ and satisfy the rit = — r/8 consistency relation. 
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Figure 11: With Na — 10^^, the figure shows the (3— Ha plane of parameter space for the AdS warp 
factor without orbifolding. The red region fits within 2a, and the yellow band is the constraint 
from COBE normahzation. We find that we need Na roughly 10^'* to fit the WMAP data. The 
blue lines are contours of r = 0.01,0.1,0.15. The model yields significantly large tensor to scalar 
ratio. The green lines are contours of 7 = 1.001, 1.5. All the models within the red region deviate 
from the slow-roll scenario, but 7 at the pivot scale is not large enough to be in the ultra-relativistic 
regime. They all belong to the intermediate region. The tightest constraint here is actually the 
bulk volume bound (4.1), with Na ^ 10^^ : no points on this figure survive the bound. 
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Figure 12: A similar graph as Figure |Tl| except that the volume of is decreased to t; = 10 
We find models that fit WMAP data, and they come with much smaller iV^ values than the case 
with V ~ 1. The graph shows a plane with Na = 20. This figure is still in the intermediate region. 



However, even if A^i is decreased significantly, no model survives the bulk volume bound (4.1) 
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Figure 13: We now take the limit Na = 1 and see if there are any models that fit the data and 
survive the bulk volume bound (4.1) at the same time. The blue lines are contours for different 
(f'pivot values. (Blue-tilt at small (ppivot-) In the 2a region of the rig constraint, i.e. the red region, 
we see that (j)pivot > 10, badly violating the bound on field range. The diamond points do show 
some points that survive the bound but they do not fit well. The figure shows that the value of 
'Ppivot decreases in the upper-right and lower-left corner, so for (3 > 10^ and f3 < 10"^, the chance 
to get small (ppivot is much larger, and these two regions are the slow-roll region and the DBl region. 
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A. The numerical approach to solving the equation of motion 



In DBI inflation, the equations of motion are solved based on the Hamilton-Jacobi ap- 
proach. The equations we want to solve are pi 



7(0) 



1 



3M2 



[y((/.) + r(7(</.)-i)] 



1 + AMjT-'^H'i^y 



(7 + 1) 



(A.l) 
(A.2) 
(A.3) 



Note that the inflaton equation of motion (^]^) (including the 4> term) is treated exactly 
here. In the Hamilton-Jacobi approach, we treat all the variables H, 7, as a function of 
(j), thus we can replace a second order differential equation ( |3.2D with a group of first order 
differential equations. The advantage of Hamilton-Jacobi approach is even if the potential 
is dependent on 7, which in turn depends on <j), one can verify that the approach still leads 
to the same equation of motion as derived from variation of the action. 
The essential idea of the numerical approach is the following. 



1. specify the initial condition at a certain position (pi; here we choose to start at the 
edge of the throat 4>r = \/TzR. We also need to give the initial value of We 
always assume that at the start of inflation, the speed of the inflaton is much smaller 
than the warped speed of light, (jii <^ T^h{(l)i)'^ and 7(0r) ~ 1. 



2. calculate H'{4)r) using ([AI2D 



^ (7^-l)T 
4M4 



(A.4) 



From ( |3.4| ), we know H' > when the brane is moving into the throat and H' < 
the brane moves out. 



3. calculate H{(Pr) using (|A.lD 



^(0) 



' [y(^) + T(7-l)] 
3M2 



(A.5) 



4. using H'{(pji) and H{(pji), calculate H{(j)ji 
method. 



numerically by the Runga-Kutta 



5. using the new Hubble parameter at — Acj)), ( |A.l ) can give "f{4>R — A(/)). If we use 
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(|A.3| ) and (|A.1| ) together, we can solve 7 through the equation 



+ 57 + C = (A.6) 



^0 



A = ^^-T (A.7) 
B = ^^^-^mV^ + 3//2-Vb + r (A.8) 



4iVA 


04 


Vo ( 




2iVA 




V^o < 






04 



C = 71^^ (A.9) 



6. use the new calculated 7(0_r — A0) in step 1 and repeat the whole process to evolve 
the equation of motion. 

When the equation of motion is solved numerically, we can get the numerical depen- 
dence of H and 7 on 0, we can further perform numerical differentiation to get H'((f)) and 
H"{cj)) and calculate all the inflationary parameters e(0), ri{(j)) and as defined in (|3.S|). 

To compare with observational data, we need to calculate the spectrum of density 
perturbation at observable scales. The scalar mode spectrum is calculated using ( 3.21| ) 



and the tensor mode using (|3.23| ). It's important to note that, because of the speed of 
sound 1/7 for the scalar mode, the tensor mode and scalar mode generated at the same (/> 
during inflation correspond to different scales when they exit the horizon. For scalar mode, 
k = aH^, while for tensor mode k = aH. In DBI scenario, generically 7 > 1, the two 
scales do not match. In order to calculate the tensor to scalar ratio correctly and compare 
with experiment, special attention must be paid to calculate the ratio on the same scale k, 
not the same inflaton value (p. To calculate the corresponding spectral index and rit, we 
avoid using any analytic formulas but simply apply numerical differentiation to the power 
spectrum. 



B. Tying the initial power spectrum to observable scales today 

The first constraint from data is the overall normalization of the scalar density perturbation. 
COBE gives us normalization at the comoving scale of the horizon today kn 

kn ~ 1.17 X IQ-^Mpc'^ 
6{kH) = (1.7 ±0.3) X 10"^ 



This value has a tilt dependence 6'^{k) = 5'^{kH) x {k/kn)^"^^ and is also calculated 
assuming no tensor contribution. 

WMAP reports its results using As{kpi^ot) = PR{kpivot) where kpi^ot = 0.002/Mpc. For 
our analysis we compare the model to the likelihood space allowed by WMAP 3 year CMB 
data [P, ^ 10, 11 1 in combination with the normalization marginalized Sloan Digital Sky 



Survey LRG matter power spectrum |12] and SNLS Supernovae sample |13]. We compare 
the DBI models to constraints from observations without imposing the consistency relation, 
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nt{kpivot) = —f{kpivot)/'^i but imposing a prior on nj, —0.02 <nt < 0., for which 

As{kpi,ot) = (23.4+;;j[^:^;) x 10-^° (b.i) 

ns{kpivot) = 0-9'''0-o.oi9(o.037) (-^-2) 

at the 68% (and 95%) c.l. and r{kpi^ot) < 0.29 at the 95% c.L. 

At the same scale the COBE and power spectrum normahzations are related by 

25 2 

■^sikpi^ot) — (kpivot) (-^-3) 

To convert kpivot today to the scale /cq at the end of inflation, we need to know the scale 
factor at the end of inflation flj^j;, a,n/;A;o = ^piuot- We can compare the energy density in 
radiation today to energy density at reheating in order to obtain ainfi- 
The radiation density is related to the photon temperature via 

Prad{a)[ML-^] = "^gnia) (^MMj%^r(a) (B.4) 

with gH{ainfi) = 106.75 compared to gniatoday) = 3.36. 

Assuming wholly efficient conversion of inflaton energy into radiation, the radiation 
density at the end of inflation is given by p'^JJ^ = H{ainfi)^- Equally today, we have 

^l,^ = ^'^ + N,nl (B.5) 
= (l + iV.((l + /.ri-l))J^° 

= (1 + 0.68) X 2.3812 x IQ-^/i-^ (T^^!^^ 

= 4.1812 x 10~^/i-2 
Hq = 8.7578 X 10~6^/iMp (B.6) 
^ladHl = 3.2069 X 10-^2° (B.7) 



If we compare this to H{ainfY from the numerical calculation then 

infl I \ / \ 4 

Prad _ 9H[ainfl) ( atoday ^ 



^rad^o gn (today) \ ai^ji 



(B.8) 



106.75 y/"^ / 3.2069 x lO^^^ox i/4 



3.36 J \ H{ainfi) 
3.17708 X 10-3° 



(B.9) 



Once we have a^n/i, we are able to calculate the scale k^ that will become the pivot scale 
kpivot today. All the observables at the pivot scale are generated at ^p^ot during inflation. 
To calculate kpivot-, we use the horizon crossing condition 

Csko = a{4>pivot)H{<j)pivot) (B.IO) 
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Combining (|B.S| ) and ( B.IO ), we get the minimum number of e- folds to solve the horizon 
problem as 

where Hj is the inflation scale and Tjih is the reheating temperature, and we assume 100% 
efficiency in converting the inflaton energy into radiation. 

We numerically interpolate at (ppivot and calculate the magnitude of scalar density 
perturbation there. For models which can obtain the normalization requirements for the 
WMAP+ SDSS LRG+ SNLS SNla constraint (Q at the 95% confidence level, we con- 
sider the behavior of the scalar spectral index, n^, tensor spectral index, nt, and tensor to 
scalar ratio, r, at (ppivot- 
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